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At a composition far above the percolation threshold, the resistance of a composite sample in- 
creases with time due to Joule heating as a constant current of sufficiently large value is passed 
through the sample. If the current is less than a certain breakdown current (/(,) the resistance 
eventually reaches a steady value with a characteristic relaxation time r^. The latter diverges with 
current / as ~ (1 — I 2 /lb 2 ) z ■ The value of the exponent z displays large fluctuations leading 
to unusual scaling of the relaxation time. It is shown that the results lead to important conclusions 
about the nature of breakdown phenomena. 

PACS numbers: 64.60.Ht, 62.20.Mk, 72.80.Tm 
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Breakdowns or fractures in random systems form 
an important class of non-thermodynamic _phase 
transitions Q. Much of the effort s HUH A 111113110 
to understand such irreversible phenomena naturally use 
the framework already developed for thermodynamic 
critical phenomena although descriptions remain far from 
complete. Consider, for example, the order of a break- 
down transition. Zapperi et al.|5j have suggested that the 
latter is first-order as a function of the external field. The 
relevant 'order parameter', elastic constant or conductiv- 
ity suffers a discontinuous change from a finite value to 
zero at a breakdown point. Using a similar but some- 
what different model, Andersen et al.Q supported the 
idea of a first order transition but only at small disorder. 
They predicted a change to a second-order transition at 
higher disorder, thus indicating presence of a tricritical 
point. On the experimental side, some recent results on 
electrical failure in composites suggest rather conflicting 
scenerios. When a sufficiently large current / is passed 
through a composite sample (a random binary mixture of 
a conductor and an insulator |l0j) with the conducting 
fraction p typically far above the percolation threshold 
p c <C p, Joule heating causes the sample resistance to 
increase. If the current is more than the breakdown cur- 
rent lb, heating eventually leads to an irreversible electric 
breakdown. It was found 8] that the resistance of such a 
sample increases as a power of time, R ~ (t^(7) — t) 
when a constant current I > i& is applied. The relaxation 
time (i.e. time-to-failure) r/j also exhibits a critical be- 
haviour (see below). However, measurements starting 
with small currents show that the sample conductance 
drops to zero as soon as the current equals or exceeds 
the breakdown current, thereby signalling a first-order 
transition. 

In order to resolve this and other issues it would be 
useful to obtain information about other relevant quan- 
tities near the transition. In this paper, we report the 
results of our dynamic measurements of thermal relax- 
ation using currents less than, but upto, the breakdown 



cent of the well-known phenomena of slowing down|ll) 
near critical points in thermodynamic systems. The re- 
laxation times also exhibits an unusual scaling relation 
due to a strong dependence on disorder. The present re- 
sults together with the earlier ones[8j constitute, to our 
knowledge, the first full description of a dynamic break- 
down phenomenon both above and below the critical pa- 
rameter. We examine below various features including 
the interplay of disorder and breakdown dynamics, which 
becomes significant particularly in the regime of weak 
disorder p c Cp<l. Few related experimental works in- 
clude investigation of avalanche dynamics |l2| and strain- 
relaxation measurements in metal networks near the per- 
colation threshold^]. 

Measurements were performed in composites of carbon 
and wax (C-W) under constant dc current condition at 
room temperature and the sample resistance was moni- 
tored as a function of time during heating. The prepa- 
ration and characterisation of the system have been de- 
scribed earlier |l4|. p c in C-W system is 0.76% (by vol- 
ume). In the tunneling regime [l5j close to p c , the resis- 
tance decreases with bias. It has been recently found 
out that there exists a Joule regime in, or a conducting 
fraction pj above, which the resistance always increases 
with bias|9j. The nominal carbon fractions of the samples 
ranging from 4.5 to 10% were above or near pj which is 
approximately 4.5%|l6j. Samples prepared with larger 
carbon fractions were mechanically unstable. On the 
other hand, samples with lower carbon fractions had tun- 
neling effect offsetting the effect of Joule heating (see Fig. 
1 of Ref. i). The geometry of current flow in a cylindcr- 
ical sample is illustrated in the inset in Fig. 1. Various 
properties of the samples used in the present work are 
given in Table I. In a typical measurement, the current 
corresponding to a certain value I would be turned on at 
time t = and the bias across the sample be measured 
at an interval of time (usually 250 msec) until it reached 
a steady value when a balance between dissipation and 
generation of heat within the sample is established. The 
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TABLE I: Various properties of the samples used, p is the 
carbon fraction by volume (%). R , p and Tho are the re- 
sistance, resistivity and relaxation time at zero current. It is 
the breakdown current and z is the exponent in Eq. Q. 



Sample 


P a 


Height 


Ro 


Po 




h 


z 


No. 




(mm) 


(O) 


(ficm) 


(sec) 


(mA) 




1 


7.5 


6.0 


422 


567 


495 


31.5 


0.175 


2 


4.5 


4.1 


287 


561 


145 


19.2 


0.187 


3 


7.5 


4.1 


139 


280 


168 


73 


0.30 


4 


4.5 


2.7 


77 


214 


455 


83.1 


0.076 


5 


7.5 


2.8 


73 


191 


113 


112 


0.18 


6 


7.5 


4.2 


04 


188 


156 


95 


0.15 


7 


7.5 


1 0x4x2' 


184 


15.2 


67 


18.7 


0.23 


8 


10 


8.3 


5.56 


5.26 


194 


135 


0.173 



"Nominal value. 

'This sample is ribbon-shaped but all other samples are cylinders 
of 10mm diameter. 
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current, not exceeding the breakdown current. The re- 
laxation during cooling from a hot steady state was also 
measured. It exhibited an anamalous behaviour. This 
will be reported elsewhere. 

Figure 1 shows a typical relaxation curve (a) during 
heating, obtained from a cylinderical sample (#3). The 
curve is well described by a simple exponential process in- 
dicated by the solid line in Fig. 1. It may be recalled [l^ 
that the relaxation function f(t) is almost universally 
given by fit) ~ exp[— (t/r) a ]. Here t is time, r is a 
relaxation time constant and a is an exponent. The 
exponential relaxation corresponding to a = 1 occurs 
mostly in simple systems such as homogeneous ordered 
solids 1 18j. A simple exponential in the present case may 
be attributed to a long-time relaxation behaviour. In 
systems with characteristic time scales such as percolat- 
ing networks, the relaxation function may change from 
a stretched exponential at short times to a simple expo- 
nential at long times [l3l Hij . For samples with large 
p ^> p c such as the ones used here, the cross-over times 
should be quite small. Close to p c , a was determined 
to be about 0.8^^ in elastic relaxation in dilute metal 
foils, and about 0.4 20] in voltage relaxation in the same 
system as the present one. 

The relaxation time in several samples of different com- 
positions and dimensions (see Table I) against current are 
shown in Fig. 2. Th of a sample was found to diverge as 
a power-law 



FIG. 1: A typical resistance relaxation curve (a) of a cylin- 
derical sample (#3) heated by passing a constant current 
of 68.5mA. Many data points have been omitted for clar- 
ity. The solid line is a fit to a simple exponential p/p = 
1.397 — 0.397 exp(—t/rh). The same data are also shown in 
log- linear plot (b). The slope of the straight line yields nege- 
tive inverse of the relaxation time Tfc=348s. The inset is an 
assembly diagram showing copper foils (c) used as electrodes 
and attached to the sample (s) with silver paints. 
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FIG. 2: Relaxation time vs. current in various samples as in 
Table I. Some data have been translated vertically by factors 
as marked. 
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Here, Tho is a prefactor, e = {I 2 /lb 2 — 1) < and z is 
a dynamical exponent. r/j , If, and z were treated as fit- 
ting parameters and its values are listed in Table I. Tho 
is really the relaxation time at zero current. It depends 
on sample dimensions (see Eq. 01 but more importantly, 
monotonically increases with disorder. In fact, r^ di- 
verges as p approaches p c . Such disorder- induced diver- 
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|9|- i? is the resistance at zero 



with R as 1^ 
current. 

The exponent z thus determined was not a constant 
even within experimental unceratinities of 5%. It ex- 
hibits an unusual sample-to-sample fluctuation as seen 
in Fig. 3 where it is plotted against a quantity r as a 
measure of disorder. (It is not practical to use p due to 
uncertainities in its values), r is defined as the ratio of 



FIG. 3: The exponent z vs. the disorder parameter r = 
Po/Pcarbon for / < /(, . The point P(l,l) is a theoretical one 
corresponding to the pure conductor. See text for furthur 
discussion. 



disorder (p — 1) and increases monotonically with disor- 
der as p decreases from 1. Within the range of disorder 
(two orders of magnitude) investigated in this work, the 
exponent has an average value of 0.184 ±0.06 i.e 33% un- 
certainities. The highest value (0.3) is about four times 
the lowest one (0.076). Without these two extreme val- 
ues, uncertainities reduce to 17% which is still beyond the 
experimental error. Clearly, the relaxation exponent z for 
I < lb is disorder-dependent. In fact, one can classify the 
breakdown-related exponents, both theoretical and ex- 
perimental, into two groups: 1) disorder-dependent and 
2) disorder- independent or 'universal'. Examples of the 
first group include, besides z of this work, the ones as- 
sociated with broken bonds in lattice models 0] while 
those in the second group include the size exponents in 
Table I of Ref. 0, roughness exponent (2lJ, relaxation ex- 
ponent z for I > It (see below). A disorder-dependent 
exponent is usually derived from measurements involv- 
ing a single sample and hence, fluctuations in its values 
possibly reflect lack of self-averaging property. 

To show how well Eq. ^ represents the relaxation time 
data, one notes that plotting (r/j/r^o) 1 ^ against I 2 /i& 2 
should lead to data collapse. This is indeed seen in Fig. 
4(open symbols, I 2 jl 2 < 1). The solid line is a fit to 
Q), which is excellent even at low currents. Thus, Eq. 
n holds not only near 1^ but also for the entire range 
of current I < If,. Goodness of the data collapse near 
lb is highlighted in the plot in the inset. Staight lines 
have slopes of unity. To complete the description of the 
dynamic breakdown, time-to-failure data (closed sym- 
bols) obtained above the breakdown point (e > 0) by 
Lamaignere et al.@ are also shown in the figure. In this 
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FIG. 4: Scaling behaviour of relaxation time below and above 
the breakdown point at I = lb. The solid and dashed lines 
are fits to (r; l /r; lo ) 1/ ' z = (=Fe) _1 respectively. The data corre- 
sponding to solid symbols are from Ref. 0. 



time at zero current as in the case e < 0. The dashed line 
is again a fit to with z = 1. The single function fits 
reasonably well all the data whereas previously, different 
functions were used to fit different ranges of data. The 
function Th ~ (I /lb — l) -2 as considered in Ref. |8jnear h 
is also shown in Fig. 4(dotted line). As seen, it does not 
fit the whole set of data. It is important to note that the 
value of the exponent z (i.e. 1) for e > is independent 
of disorder in contrast to the varying exponent for e < 0. 
This is ensured by Eq. ^ and the requirement that Th 
must vary as I~ 2 at large currents || as dissipation of 
heat becomes negligible at large currents. 

It is now useful to summarize characteristics below and 
above the breakdown point / = lb- (i) For e < 0, the re- 
sistance ratio R/R (see, for example, Fig. 1 of Ref. 0) 
smoothly increases with current to a limiting value (or 
breakdown ratio) T « 1.5[l6J at e = 0. R = R(t — > oo) is 
the steady resistance at a given current. It is best fitted 
by an expression R/R = T + a\e + a-2e 2 + a^e 3 that has 
a finite slope a\ at e = 0. For e > 0, there is no steady 
state. On the other hand, 77, diverges as a power-law of e 
on both sides; ii) z is not same on both sides as discussed 
above. This is in contrast to equality of the exponents 
on both sides, that is expected from scaling hypotheses 
in both static and dynamic critical phenomena in ther- 
modynamic systems; iii) While the scaling in Fig. 4 is 
itself an expression of 'universality' 0, the fluctuations 
in z for e < renders the notion of an universality class 
unteneble; iv) At e < 0, both R and Th are reversible 
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versible change in the micro-structure should also lead to 
irreversibility in those quantities. This is incompatible 
with the picture that emerges in breakdown models with 
annealed disorder^ [3- 0- 0- where the global break- 
down is preceded by increasing bursts of irreversible bond 
breaking. The composite samples can be considered be- 
longing to a system with quenched disorder for e < 0, 
and a system with annealed disorder for e > 0. In view 
of all these, it is difficult to avoid the conclusion that the 
usual thermodynamic classification of phase transition is 
inadequate to describe breakdown phenomena. 

The presence of disorder is expected to have a profound 
effect on the breakdown processes. One of its manifesta- 
tion lies in the conceptual difficulties of taking the limit 
of disorder going to zeroQ. This is aptly illustrated in 
the present case. Let us first consider divergence of the 
relaxation time in a homogeneous medium corresponding 
to zero disorder. In fact, it has been long discussed [22} 
albeit couched in the language of stability. The Joule 
heating under a constant current / in a medium with a 
positive temperature coefficient of resistance (3 leads to 
either a steady state or breakdown depending on whether 
the heat generated is removed quickly enough or not. Full 
solutions of heat flow have been obtained in some regular 
geometries [22} where is given by 
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l 2 {Ka-bI 2 T~l 2 K{\-I 2 /I b 2 ) 
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where lb — na/b, k is thermal diffusivity, and I is the 
smallest distance for flow of heat, a, b are constants that 
depends on the boundary conditions and various mate- 
rial constants including (3. Clearly, z = 1 in the case of 
zero disorder (p — 1). In passing let us note that ac- 
cording to Eq. I 2 rather than I is the proper variable 
to use in problems involving the Joule heating. Consid- 
eration of temperature-coupled resistance is an essential 
ingradient in obtaing (J2J ■ Its absence in the dynamic fuse 
model of Sornette et al.@| gives rise to a relaxation time 
independent of current. 

As the disorder is reduced (p — > 1), the relaxation ex- 
ponent increases from ~ 0.18 to 1. Thus, the dependence 
of z on disorder is much stronger than the fluctuations 
in its values suggest. Presently, there is no theory of z 
for p < 1. An interesting question arises at this point 
as to how z from a lower value at higher disorder would 
approach 1 at the point P in Fig. 3. Three possible 
curves are shown schematically. The curve a is intuitive 
and has a zero slope at P. This means that the sys- 
tem having a small disorder can be simply considered as 
a homogeneous one with an effective thermal diffusivity. 
This may be possible if the melting point of the conduc- 
tor is less than that of the insulating matrix. The curve 
c corresponds to the situation where even a slight dis- 
order leads to an abrupt fall in z. This may happen if 



at P. Clearly, furthur theoretical efforts are necessary for 
understanding of z as a function of disorder. Similar 
problem^ concerning the breakdown ratio T exists near 
p~l. " 

In conclusion, a detailed description of the dynamic 
electric breakdown driven by external current, both be- 
low and above the breakdown current, was presented. 
It was shown that the classification of thermodynamic 
phase transitions is inadequate for breakdown transi- 
tions. The present models having annealed disorder lack 
reversibility found in composites. 

We acknowledge several discussions with Bikas 
Chakrabarti and the assistance of Arindam Chakrabarti 
in the preparation of C-W samples and processing data. 
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